Rapid identification and interpretation of gene–environment associations using the new R.SamBada landscape genomics pipeline

sam βada is a genome–environment association software, designed to search for signatures of local adaptation. However, pre‐ and postprocessing of data can be labour‐intensive, preventing wider uptake of the method. We have now developed R.SamBada, an r ‐package providing a pipeline for landscape genomic analysis based on sam βada , spanning from the retrieval of environmental conditions at sampling locations to gene annotation using the Ensembl genome browser. As a result, R.SamBada standardizes the landscape genomics pipeline and eases the search for candidate genes of local adaptation, enhancing reproducibility of landscape genomic studies. The efficiency and power of the pipeline is illustrated using two examples: sheep populations from Morocco with no evident population structure and Lidia cattle from Spain displaying population substructuring. In both cases, R.SamBada enabled rapid identification and interpretation of candidate genes, which are further discussed in the light of local adaptation. The package is available in the r CRAN package repository and on GitHub (github.com/SolangeD/R.SamBada).     

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species‐wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well‐watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA‐seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species‐wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population‐specific responses. Weighted gene co‐expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Patterns of transposable element variation and clinality in Drosophila

Natural populations often exist in spatially diverse environments and may experience variation in the strength and targets of natural selection across their ranges. Drosophila provides an excellent opportunity to study the effects of spatially varying selection in natural populations, as both Drosophila melanogaster and Drosophila simulans live across a wide range of environments in North America. Here, we characterize patterns of variation in transposable elements (TEs) from six populations of D. melanogaster and nine populations of D. simulans sampled from multiple latitudes across North America. We find a nearly twofold excess of TEs in D. melanogaster relative to D. simulans, with this difference largely driven by TEs segregating at the lowest and highest allele frequencies. We find no effect of latitude on either total TE abundance or average TE allele frequencies in either species. Moreover, we show that, as a class of mutations, the most common patterns of TE variation do not coincide with the sampled latitudinal gradient, nor are they consistent with local adaptation acting on environmental differences found in the most extreme latitudes. We also do not find a cline in ancestry for North American D. melanogaster—for either TEs or single nucleotide polymorphisms—suggesting a limited role for demography in shaping patterns of TE variation. Though we find little evidence for widespread clinality among TEs in Drosophila, this does not necessarily imply a limited role for TEs in adaptation. We discuss the need for improved models of adaptation to large‐scale environmental heterogeneity, and how these might be applied to TEs.     

Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes

Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

基于SPEI指数的淮河流域干旱时空演变特征及影响研究

运用淮河流域149个气象站1962—2016年逐日气温、降水资料以及历史旱情资料,基于SPEI、EOF和M-K等方法分析淮河流域的干旱时空特征,研究干旱的时空演变规律并揭示其对农业生产的影响。结果表明:(1)基于SPEI得到的干旱频次与受灾、成灾面积的相关性通过了0.1的显著性水平检验,表明SPEI在淮河流域具有较好的适用性;(2)淮河流域干旱发生时间差异明显,干旱次数呈现波动变化,发生重旱和特旱次数占总干旱的比重是20.0%,其中重旱和特旱在1960s比重最大(24.8%),其次是2010s(15.8%),在1980s比重最低(10.0%);(3)干旱的空间分布差异大,淮河流域干旱频率在27.76%—36.04%之间,西北部和东南部发生干旱强度较西南部、东北部及中部低;(4)淮河流域总体呈干旱化的趋势,从中部到四周呈现由高到低递减的趋势变化,且空间模态表现为全区一致型、南北相反型和东西相反型。     

贵州喀斯特区C4植物物种组成与水分生态型划分

该研究通过查阅文献、核对贵州大学林学院标本库及现场群落调查与标本采集,并运用碳同位素比值法研究典型C_4植物水分利用特性,探索C_4植物在喀斯特植被恢复中的地位,进而揭示贵州喀斯特地区C_4植物资源的基本特征。结果表明:贵州喀斯特区共有C_4植物141种,隶属于74属15科,分别占全国科属种的62.50%、46.25%、24.48%,以禾本科(Gramineae)和莎草科(Cyperaceae)为主;区内C_4植物种均为一年生或多年生草本,多年生植物种略多于一年生植物种;水分生态型整体偏旱生,旱生和中生植物分别占总数的24.82%和31.21%;喀斯特区C_4植物具有高水分利用效率,但不同水分生态型间差异不显著;贵州喀斯特区C_4植物资源具有资源丰富、利用途径广泛、能长期利用、竞争力强、能大面积分布、偏旱生且水分利用幅度广的基本特征,适合喀斯特区生境,自然状态下多为恢复早期物种,有利于喀斯特区生态恢复。在贵州喀斯特恶劣生态环境下C_4植物有较好的生态适应性,并表现出较高的药用、食用、饲用、景观应用等价值,对其开发利用对贵州经济、社会发展及生态恢复有重要意义。     

海州湾鱼类群落平均营养级和大型鱼类指数的变化特征

为探究近年来海州湾鱼类群落的变化特征,本研究根据2011年及2013-2017年春季(5月)和秋季(9-10月)在海州湾及其邻近海域进行的渔业资源底拖网调查,分析了海州湾鱼类群落平均营养级(MTL)和大型鱼类指数(LFI),对海州湾鱼类群落结构特征进行研究.结果表明:海州湾的优势鱼种主要有大泷六线鱼、方氏云鳚、尖海龙、小黄鱼、长蛇鲻等,且优势鱼种季节性变化明显.海州湾鱼类群落的平均营养级存在明显的年际和季节变化,总体上秋季的MTL高于春季,而且秋季MTL变化具有滞后性.LFI计算结果表明,近年来海州湾及其邻近海域大个体鱼类资源量有所减少,鱼类群落结构呈现出明显的小型化趋势.